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i We provide a detailed study of the phenomenology of orbifold compactifications of the het- 

' erotic string within the context of supergravity effective theories. Our investigation focuses on 

those models where the soft Lagrangian is dominated by loop contributions to the various soft 
supersymmetry breaking parameters. Such models typically predict non-universal soft masses 
and are thus significantly different from minimal supergravity and other universal models. We 
consider the pattern of masses that are governed by these soft terms and investigate the im- 
plications of certain indirect constraints on supersymmetric models, such as flavor-changing 
neutral currents, the anomalous magnetic moment of the muon and the density of thermal relic 
neutralinos. These string-motivated models show novel behavior that interpolates between the 
phenomenology of unified supergravity models and models dominated by the superconformal 
anomaly. 
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Introduction 



The recent interest in models where supersymmetry breaking is transmitted from a hidden sector to 
our observable world via the superconformal anomaly [1, 2, 3] has served as a reminder that there 
are occasions where working at one loop order is more than a theoretical luxury, but an absolute 
necessity if one is to understand even the gross features of some models. Recently the complete 
one loop supergravity correction to soft supersymmetry breaking terms has been obtained [4, 5], 
a subset of which arc the aforementioned "anomaly-mediated"' terms. The broad features of these 
corrections, in particular for gaugino masses, was studied in a heterotic string context in [6]. In the 
current work we wish to study in greater detail the low energy phenomenology of models based on 
orbifold compactifications of the weakly-coupled heterotic string in which these loop corrections are 
significant. To that extent wc will choose two "laboratories" in which tree level soft supersymmetry 
breaking terms are either absent or greatly suppressed. 

Our first example will involve an observable sector in which the matter fields are assumed to 
be from the untwisted sector of the orbifold compactification. In the event that supersymmetry 
breaking is transmitted by the compactification moduli T", whose vacuum expectation values 
determine the size of the compact manifold, then this model will exhibit a no-scale pattern of soft 
terms with vanishing soft supersymmetry breaking terms at tree level. Clearly the phenomenology 
of such models will depend critically on the form of the loop induced soft terms. Indeed, the 
"sequestered sector" models considered in [2] were of just such a form. 

In our second example we will switch our focus to cases in which it is the dilaton field S, whose 
vacuum expectation value determines the magnitude of the (unified) coupling constant Qstr at the 
string scale, that participates exclusively in supersymmetry breaking. We will work in the context 
of models in which string nonperturbative corrections to the Kahler potential act to stabilize the 
dilaton in the presence of gaugino condensation [7, 8]. Such theories predict a suppression of dilaton 
contributions to soft terms relative to those of the supergravity auxiliary field - thereby suppressing 
tree level contributions to supersymmetry breaking relative to certain loop contributions. 

After describing our framework and introducing both the tree level and one loop soft supersym- 
metry breaking terms in Section 1, we introduce the parameter space for our two broad classes of 
models in Section 2. The subsequent section describes the analysis tools we employ and the types of 
low-energy phenomena we will study before describing the implications of low-energy experimental 
constraints on each of our models in Section 4. 

In general we find that these examples have a phenomenology that combines the features of 
minimal supergravity (mSUGRA) models with those of anomaly-mediated models. Parameters 
related to the orbifold compactification and the stabilization of string moduli interpolate between 
these regimes. Current observations and particle mass limits from collider experiments already shed 
light on the nature of moduli stabilization. Future discoveries of superpartners will likely provide 
strong evidence in favor of one stabilization mechanism or another in the context of weakly-coupled 
heterotic string models. 



2 



1 The structure of heterotic orbifold models at one loop 



We work in a framework in which the chiral superfields Z can be divided into two classes: 
observable sector superfields, noted Z*, charged under the observable sector gauge symmetries; 
and hidden sector superfields, noted Z". Since our interest here is primarily on broad issues of 
phenomenology, the hidden sector fields we will consider are the dilaton S in the chiral multiplet 
formulation and the three diagonal Kahler moduli T". 

In the orbifold compactifications wc study here, the tree-level Kahler potential for the moduli 
and the matter superfields is known. For the moduli sector we have 

3 

if(5,5;r",T") = + - ^ln(r"+T") , (1.1) 

a=l 

where we prefer to leave the form of the dilaton Kahler potential k unspecified at this point. Its 
precise form depends on how one stabilizes the dilaton. We will return to this issue when we discuss 
the dilaton-dominated scenario of Section 2.2. As for the observable sector matter fields with 
modular weights raf associated with each of the three T", we assume a diagonal Kahler metric 
given by 

K,- = K,[Z-)b,j + 0{\Z'W (1.2) 

with 

,ci(z") = [](r" + T")"?. (1.3) 

a 

In the interests of simplicity, we assume that the three Kahler moduli can be treated as equivalent 
so that 

i^(5,5;T,T) = A:(5 + 5)-31n(r + T); Ki = (T + T)"S (1.4) 

where n-i = ^„ nf. 

The tree-level gauge kinetic functions faiZ"-), one for each gauge group Qa, are given in the 
weak coupling regime by 

fain = S. (1.5) 

Their vacuum expectation values give the associated gauge couplings (Re/a) = l/g^- To (1.5) will 
be added certain string threshold corrections when we exhibit the one loop soft supersymmetry 

breaking terms below. 

The scalar potential, written in terms of auxiliary fields, is given by the expression^ 

V = Kj-jF^F"^ - \mM, (1.6) 
^We will assume vanishing D-terms in what follows. 
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with Kj-j = d^K/dZ^dZ being the Kahler metric. The fields in (1.6) are the auxihary fields 
associated with the chiral superfields Z^ while the field M is the auxiliary field of the supergravity 
multiplet. Solving the equations of motion for these auxiliary fields yields 

= -e^/^K^^ (Wj^ + Kjfw), (1.7) 

M = -Se^/'^W, (1.8) 

with being the inverse of the Kahler metric. Note that these expressions are given in terms 

of reduced Planck mass units where we have set Mpl/x/Stt = 1. The supergravity auxiliary field is 
related to the gravitino mass by 

M3/2 = -\ (m) = (e^/W) . (1.9) 

We adopt the ansatz of Brignole et al. [9] in which one assumes that the communication of 
supersymmetry breaking from the hidden sector to the observable sector occurs through the agency 
of one of the moduli - in this case either the dilaton S or the (universal) Kahler modulus T - by 
the presence of a non- vanishing vacuum expectation value of their auxiliary fields F^ or F-^. In 
principle both types of moduli could participate in supersymmetry breaking, and so one typically 
introduces a "Goldstino angle" 9 to parameterize the degree to which one sector or the other feels 
the supersymmetry breaking. While this is a common practice, all known explicit models of moduli 
stabilization in the heterotic string context predict that this "angle" will be some integer multiple 
of tt/2; that is, only one of the two classes of moduli participates in supersymmetry breaking. 

If these are the only sectors with non-vanishing auxiliary fields in the vacuum, then the further 
requirement that the overall vacuum energy at the minimum of the potential (1.6) be zero allows 
us to immediately identify (up to phases, which we set to zero in what follows) 

F^ = — ]=Mk;^^^^ sine = V3M3/2k;^^^^ sine, (1.10) 
v3 

= — ^MK~^^^'^ cose = V3M3/2K~^^^'^ cose, (1.11) 
v3 

where the last equality holds for the vacuum expectation values, using (1.9) above. We should note 
that the condition of vanishing vacuum energy is a necessary one to employ the above parameteri- 
zation. In this work we will avoid discussing specific models of dynamical supersymmetry breaking 
and moduli stabilization, but any such model must include some mechanism for engineering a 
vanishing vacuum energy in order to make contact with the results presented here. 

1.1 Modular invariance and tree level supersymmetry breaking terms 

The soft supersymmetry breaking terms in string-derived supergravity depend on the moduli 
through the observable sector superpotential and this is, in turn, determined by insisting on mod- 
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ular invariance of the low-energy effective Lagrangian. The diagonal modular transformations 

„ aT — ib 



icT + d' 



ad-bc=l, a,b,c,deZ, (1-12) 



leave the classical effective supergravity theory invariant, though at the quantum level these trans- 
formations are anomalous [10, 11, 12, 13]. This anomaly is cancelled in the effective theory by the 
presence of a universal Green-Schwarz counterterm and model-dependent string threshold correc- 
tions [14, 15], which we describe below. 

A matter field of modular weight nj transforms under (1.12) as 

^ (zcT + df'Z* (1.13) 

while the Kahler potential of (1.4) undergoes a Kahler transformation K ^ K + 3{F + F), with 
F = \n{icT + d), under (1.12). Therefore the classical symmetry will be preserved provided the 
superpotential transforms as [16] 

W ^ W {icT + d)~^ . (1.14) 

To ensure this transformation property the superpotential of string-derived models has a moduli 
dependence of the form 

w,^k = w,,k [^(r)]-2(3+-.+"^+"'=) . (1.15) 

where Wijk = d^W{Z^)/dZ''dZ^ dZ^ . The function ry(T) is the classical Dedekind eta function 

oo 
n=l 

and it has a well-defined transformation under (1.12) given by 

r?(T) ^ {icT + d)^/^ ri{T). (1.17) 
We will also have to introduce the modified Eisenstein function 

G2(t,t)^2C(t) + ^-, where C(T) = (1.18) 

which vanishes at the self-dual points t = 1 and t = e^'^/^. 

We are now in a position to give the tree level soft supersymmetry breaking terms. The tree 
level gaugino mass for canonically normalized gaugino fields is simply 

= ^F-dX- (1-19) 
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We define our trilinear A-terms and scalar masses for canonically normalized fields by 

= ^ E Mj^^'^l'^Wi^^.z'z^z^ + h.c. = \Y, AijkeK/2(KiACjACk)-'/'Wijkz'zjz'^ + h.c. , (1.20) 



6^ '-"^ 6 4r 

i]k ijk 

where = k^^z'^ is a normalized scalar field, and by 



VM = YlM!'^i\^f = YlMi\zf. (1.21) 

i i 

With these conventions our tree level expressions are 

= {F^dnHmKjKke-'' /Wijk)) ■ (1.22) 



(M°)2 = (^^-F-F^dndmlnK,). (1.23) 

If we specialize now to the case of (1.10) and (1.11) with modufi dependence given by (1.4), (1.5) 
and (1.15), then the tree level gaugino masses (1.19), A-terms (1.22) and scalar masses (1.23) 
become 

a 2 

= i3 + ni + n, + nk)G2{t,i)F^ -ksF^ 



MM 

'■'(tT't) 



The expressions in (1.24) are to be understood as vacuum values, though we will drop the cumber- 
some brackets (...) from here onwards. 

These tree level soft terms have been studied extensively in the past and we wish to understand 
what impact the less-studied loop corrections can have on the phenomenology of orbifold models. 
We can now appreciate the importance of the two cases we intend to study in Sections 2.1 and 2.2. 
In the case where the modular weights of the observable sector are universally = — 1, as would be 
the case for the untwisted sector of orbifold compactifications, then from (1.15) it is clear that the 
moduli do not couple to the observable sector through the superpotential and we have a true no- 
scale model. Indeed, substituting (1.11) into (1.24) we find that if the dilaton does not participate 
in supersymmetry breaking all {6 = 0) tree level soft terms are precisely zero in such a case. This 
is the situation we will investigate in Section 2.1. 

Alternatively, if the dilaton is the primary source of supersymmetry breaking in the observable 
sector, loop-level corrections to at least the trilinear A-terms and gaugino masses will continue to 
be important provided the magnitude of F^ is suppressed relative to that of M, as it is in models 
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where nonperturbative corrections to the Kahler potential are used to stabilize the dilaton [17]. 
In such a scenario the modular weights of the matter fields are irrelevant at the tree level, so we 
will continue to assume = — 1 for the sake of convenience when we investigate these models in 
Section 2.2. We thus need to turn our attention to the study of soft supersymmetry breaking at 
one loop. 



1.2 General one loop corrections to soft supersymmetry breaking 



In this section we aim to provide sufficient background to justify the form of the one loop soft term 
expressions which are our goal, as well as explain some notation wc will need for our phenomeno- 
logical analysis. For more complete explanations of what is contained in this section, including ex- 
pressions with arbitrary modular weights and three independent Kahler moduli, the reader should 
consult the precursor to this work [6]. 

Wc begin with gaugino masses which can be understood as a sum of loop-induced contributions 
from the field theory point of view, and terms that can be thought of as one loop stringy corrections. 
The field theory loop contribution can be derived completely from the superconformal anomaly and 
is given by [4, 18] 



47r2 



(1.25) 



where Cq, C* are the quadratic Casimir operators for the gauge group Qa in the adjoint represen- 
tation and in the representation of Z*, respectively, and ba is the beta-function coefficient for the 
group Qa- 

ha = -^\^Ca-yCl\. (1.26) 



167r2 



As mentioned in the previous section one expects modular anomaly cancellation to occur 
through a universal Green-Schwarz counterterm with group- independent coefficient (5gs- Such a 
term can be thought of as a loop correction that contributes to gaugino masses in the form 

M So 



Jgs 



(1.27) 



2 {t + t) 167r2 ■ 

In addition string threshold corrections generally appear in the effective theory, which may be 
interpreted as one loop corrections to the gauge kinetic functions of the form 

<^GS 



/i(Z")= In r?2(T) 



167r^ 



+ ba 



which generate loop contributions to gaugino masses given by 



Mi 



Ith 



167r2 



(1.28) 



(1.29) 
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Putting together the tree level gaugino mass with the loop contributions (1.25), (1.27) and (1.29) 
gives 



Ma 



{2 ^ + ba G2it,i)F^ + ^baM+[l-2bX]F'^^ (1.30) 
where we have defined the quantity 

To understand the form of the one-loop A-terms and scalar masses it is necessary to describe how 
field theory loops are regulated in supergravity, seen as an effective theory of strings. The regulation 
of matter and Yang-Mills loop contributions to the matter wave function renormalization requires 
the introduction of Pauli-Villars chiral superfields = $*, and which transform according 
to the chiral matter, anti-chiral matter and adjoint representations of the gauge group and have 
signatures r]A = respectively. These fields are coupled to the light fields through the 

superpotential 

W{^^, Z') = ^Wij{Z^)^'^^ + V2$"$i(r„Z)^ + • • • , (1.32) 
where is a generator of the gauge group, and their Kahler potential takes the schematic form 

i^Pv = E 4(^^)1^^!', (1-33) 
A 

where the functions ka are a priori functions of the hidden sector (moduli) fields. These regulator 
fields must be introduced in such a way as to cancel the quadratic divergences of the light field loops 
- and thus their Kahlcr potential is determined relative to that of the fields which they regulate. 

The PV mass for each supcrfield <1>'^ is generated by coupling it to another field 11^ = (11*, 11*, 11") 
in the representation of the gauge group conjugate to that of through a superpotential term 

Wm = Y.^'A{Z'')^''Jl'', (1.34) 
A 

where ij,a{Z^) can in general be a holomorphic function of the light superfields. There is no 
constraint on the Kahler potential for the fields 11"^ as there was for those of the However, if 
we demand that our regularization preserve modular invariance then we can determine the moduli 
dependence of the regulator fields 




Kf = Ki = {T + T)^^, 

4 = ga'e^ = ga^e'^{T + T) 



* - (1.35) 
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and we can determine the moduli dependence of the supcrsymmctric mass IjLa{Z^) in (1.34) 

fiAiZ"") = fiA{S) [^(r)]-2(3+"^+9^) , (1.36) 

with HA and qa being the modular weights of the fields and 11"^, respectively. Furthermore, we 
can postulate the form of the moduli dependence of ka for the mass-generating fields 

: hi^ = hAiS + S){T + T)i^. (1.37) 

So at this point the dilaton dependence in the superpotential term (1.34) and the functions Ha, as 
well as the modular weights qa of the fields 11"^, are new free parameters of the theory introduced 

at one loop as a consequence of how the theory is regulated. Given (1.34) we can extract the 
Pauli-Villars masses that appear as regulator masses in the logarithms at one loop 

ml = e^{Kt)-'/\ni)-'/'M', (1.38) 

with mA = {rrii, rhi, rua) being the masses of the regulator fields respectively. 

In terms of these regulator masses, the complete one loop correction to the trilinear A-terms 
and scalar masses in a general supergravity theory was given in [5]. In this survey we will simplify 
the expressions and reduce the parameter space by making some reasonable assumptions. Let us 
first assume that the functions ij,a{Z^) that appear in (1.34) and (1.38) are proportional to one 
overall Pauli-Villars scale fipy so that fii = = ^^i = /^pv- This scale is presumed to represent 
some fundamental scale in the underlying string theory. Let us further assume that there is no 
dilaton dependence of the PV masses so that hA{S + S) is trivial and /Xpv is constant. With these 
simplifications the complete trilinear A-term at one loop is given by 

= ^4fe-^7iM-G2(i,t)F^(^E7"Pm + E7i>;rn) 



-In [{t 



\ a Im / 

+2^7fM0ln(A/^R) +E7i™4mln(A//^R) + cyclic(ijfc), (1.39) 

a Im 

where we have defined the following combinations of modular weights from the Pauli-Villars sector 

Pij =^ + ^{ni + nj+qi + qj), Pia = ^{^ + Qa + Qi - ni) (1.40) 

which we will refer to as "regularization weights" in reference to their origin from the PV sector of 
the theory.^ 

^Note that as we are considering the case with universal modular weights rii = —1, these expressions can be 
reduced to pij = 2+ ^{qi + qj) and Pia = 2 + ^{qa + qi). 
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In (1.39) and A'^^^^^ are the tree level gaugino masses and A-terms given in (1.24) and the 
parameters 7 determine the chiral multiplet wave function renormalization 



327r2 



kl 



(1.41) 



We have implicitly made the approximation that generational mixing is unimportant and can be 
neglected in (1.39), and that motivates the definitions 



jk a. 



a ^ _9o^(rp2Y jk ^ 



-,K 



■^^(KiKjKk) \Wijk\ . 

The scalar masses are obtained similarly and take the form 



(1.42) 



{M^y 



jk 



+ h.c. 



1 



+ { F^G, {t, t) I E ^fPiaM'a + ^ E ^i PM% I + h.c. 



■In 



+ E7^V K)^ + + (4^)^ 

+ E7f [3(M°)2-(M°)2] ln(M?,/Mli) 

a 

jk 



(1.43) 



with M° being the tree level scalar masses of (1.24). 

To put these expressions into a less cumbersome form we will consider the case where the various 
regularization weights pia and pjk can be treated as one overall parameter p. Then inserting the 
tree level soft terms (1.24) into (1.39) and (1.43) yields these soft terms in the final form we will 
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use in the following analysis: 

Aijk = -^F^-^7«M-p7,G2(i,tlF^+7ii^^{ln(/x2j/x|)-pln[(t + t)|r?(i^ 

+cyclic(ijA;) (1-44) 

r |M|2 



Mt 



+(l-p)7i 



IMI' 



MF^ 



6 



+ h.c. \ + <^ p7iG2 i) 



F^F^ 



+ h.c. 



where 7^ is a shorthand notation for 



-Jlltgl + ^.^.-E^f ) (ln(A/Ml) - pin [(t + 1^(01']) 

a jk 



a jk 



(1.45) 



(1.46) 



The adoption of one overall regularization weight p makes it possible to identify the quantity 
ln(|Ltpy//x|j) — pin \{t + f)|?7(t)|^] as a stringy threshold correction to the overall PV mass scale, or 
effective cut-off, ;Upv The one loop expressions of (1.30), (1-44) and (1-45) will be our starting point 
in the following two sections where we investigate their phcnomcnological significance in scenarios 
where tree level soft supersymmetry breaking is suppressed or vanishing. 



2 Classification of heterotic orbifold models 

As a result of the assumptions made in Section 1 our general parameter space is defined by seven 
quantities, if we assume a tree-level Kahler potential for the dilaton of the form k{S, S) = — ln(6'+5') 
with {s + s) = 2/^1^5^ ~ 4. Three of these parameters represent scales: the gravitino mass M3/2, 
the boundary condition scale yUuv at which the soft supersymmetry breaking terms appear and the 
Pauli-Villars cutoff scale yUpy. The other four are parameters describing the moduli sector: the 
vacuum expectation value of the real part of the (uniform) Kahler modulus {t + f), the modular 
weights of the PV regulators parameterized by p, the value of the Green-Schwarz coefficient (5gs 
and the Goldstino angle 9 which determines the degree to which the dilaton and moduli F-terms 
participate in the transmission of supersymmetry breaking. 

2.1 Moduli dominated scenarios 

This section assumes that i) all fields have modular weight ni = —1, which gives a no-scale structure, 
and ii) we have moduli-dominated supersymmetry breaking {9 vanishes, cos ^ = 1) so that the tree 
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level soft terms are precisely zero. In addition, those loop-induced soft terms in (1.45) proportional 
to tree level quantities are also zero, leaving only those arising from the superconformal anomaly 
and those related to the Pauli-Villars sector: 



167r2 + 



G2 (t,i)F^ + ^6aM|, 



Aijk = ~jiM -jrfiG2{t,t)F'^ + cjdic{ijk), 

IMP 

Mi = (I-pH^-^. (2.47) 

Note that in the special case where the moduli are stabilized at one of their two self-dual points 
t = 1 and t = e'"^/^ we recover the universal "anomaly-mediated" results 

Ma = 5a()")^Y 

^M 

Ajk = - in + 7j + Ik) Y 

Mf = (2.48) 

The scenario that has come to be referred to as the "Anomaly-Mediated Supersymmetry Breaking" 
(AMSB) scenario [1, 2, 3], in which scalar masses first appear only at two loops, is obtained in the 
limit as the phenomenological parameter p approaches unity. This would be the case, for example, 
if the mass-generating Pauli-Villars fields 11^ had the same Kahler metric as the regulator fields 
and thus qi = —1, qi = 1 and qa = —S. In any event, it is clear that the parameter p defines a family 
of "anomaly mediated" models for the case of moduli stabilized at self-dual points. Importantly, 
for p < 1, as would be typical in orbifold models, the anomaly-mediated masses for the matter 
scalars are positive with only the Higgs fields having negative squared masses at the scale /^uv- 

The contrast between these results and those of the standard AMSB case was discussed in [5, 6], 
but it is important to point out that the vanishing of the scalar masses at one loop in the p ^ 1 limit 
clearly demands a two-loop analysis of these soft terms. A treatment of supergravity radiative cor- 
rections, regularized by the Pauli-Villars mechanism, at two loops would be a massive undertaking. 
Some subset of these corrections would undoubtedly be the soft terms found in [2, 3] using a spurion 
technique. As is now well known, these two-loop terms arising from the superconformal anomaly 
would imply negative squared masses for the scalar leptons. To remedy this situation, the AMSB 
model, as it has been institutionalized in studies such as the Snowmass Points and Slopes [19], 
adopts an otherwise ad hoc scalar mass contribution that is universal and sufficiently large. Not 
surprisingly, these models are often characterized as having light sleptons as a key feature of the 
phenomenology. 
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While several extensions of the original (or "minimal" ) AMSB model now exist which address 
the slepton mass problem in a variety of ways,^ the generalized anomaly mediated model we present 
here suffers from no such problem when p ^ 1. While the phenomenology associated with the 
gaugino mass spectrum is familiar from past studies of the AMSB paradigm, the scalar mass sector 
looks quite different. We will begin our analysis in Section 4.2 with this generalized anomaly 
mediated model (which we will refer to as the "PV-AMSB model" ) before enlarging the parameter 
space to include all moduli-dominated models in Sections 4.3 and 4.4. 



2.2 Dilaton dominated scenarios with Kahler suppression 

We next turn to the opposite extreme, where the dilaton is the primary source of supersymmetry 
breaking in the observable sector (sin 9 = 1). In such a scenario we would ordinarily expect the loop 
corrections such as those in (1.25), (1.39) and (1.43) to be small perturbations on the dominant 
(and universal) tree level contributions and thus safely neglected. The phenomenology of the 
dilaton-domination scenario at tree level has been extensively studied in the literature [21]. 

The bulk of these studies, however, have considered only the case of the standard Kahler 
potential k{S, S) = — ln(5 + 5) derived from the tree level string theory. It is well known that when 
this form is used, stabilizing the dilaton at acceptable weak-coupling values is extremely difficult. In 
the case where a dilaton potential is imagined to arise through field theory nonperturbative effects 
(such as gaugino condensation) it has been shown that one condensate alone can never succeed in 
stabilizing the dilaton in such a way as to reproduce the weak coupling observed in nature [22, 23]. 
Furthermore, while multiple condensates can be used to provide the requisite stabilization [24, 25], 
these models tend to predict a moduli-dominated scenario such as those in Section 2.1. 

In contrast, if one postulates a nonperturbative correction of stringy origin to the dilaton Kahler 
potential, as was first motivated by Shenker [26], then one condensate can indeed stabilize the 
dilaton at weak coupling while simultaneously ensuring vanishing vacuum energy at the minimum 
of the potential. The power of this approach was first demonstrated in an explicit model of modular 
invariant gaugino condensation in [7, 8] and confirmed in [27], providing a concrete realization of 
the so-called "generalized dilaton domination scenario" [23]. 

The key feature of such models is the deviation of the dilaton Kahler metric from its tree level 
value. If we imagine the superpotential for the dilaton having the form W{S) oc 6"^"^/^^+, with 
being the largest beta-function coefficient (1.26) among the condensing gauge groups of the hidden 
sector, then it is clear from (1.7) and (1.8) that requiring the potential (1.6) to vanish implies 



k 

264 



3 ^ {k'')-^/'^ = VS ^ — . (2.49) 

1 - o6_l/c,. 



{k' 
is i 

merely of requiring a vanishing vacuum energy in the dilaton-dominated limit. The explicit model 



The condition in (2.49) is independent of the means by which the dilaton is stabilized and is a result 



*For some early solutions, see [20]. 
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of [8], however, was able to achieve precisely this relation with {kg) = —g1,^^/2, using a correction 
to the dilaton action that involved tuning 0(1) numbers only.^ 

We can parameterize the departure that (2.49) represents from the tree level Kahler metric 
(^s?^)^^^) = (l/(* + ^)) = 9str/^ — 1/4 by introducing the phenomenological parameter 



1/2 



jLtrue 



so that the auxiliary field of the dilaton chiral supcrmultiplct can be expressed as 



(2.50) 



(2.51) 



The importance of (2.49) for our purposes is to recognize that the factor of 6_|_, containing as it 
docs a loop factor, will suppress the magnitude of the auxiliary field _F'^ relative to that of the 
supcrgravity auxiliary field M through the relation (2.51). Therefore tree level soft terms for 
the gaugino masses and trilinear A-terms will be of comparable magnitude to the loop-induced 
soft terms arising from the superconformal anomaly (though scalar masses will be only negligibly 
altered from their tree level values). We are thus led to consider the phenomenology of models 
given by the following pattern of soft supersymmetry breaking terms: 

5a (m) 



Ma = 



A. 



ijk 



Mf = 



2 

|M|2 



haM+[l-2b'aks] 

- ^7iM + jiF^ {ln{p,ljf^%) -pin [{t + f)|r/(t)|^] } + cyclic(ij/c) 

(Hi,lJi,j,)-pln [{t + tMt)\']) 



l + 7.-|E'^-2E^' 

- jk 



\^ MF^ 



(2.52) 



where we have dropped terms of O (l/(167r^)^) in the scalar masses. These models will be studied 
in Section 4.5. 



3 Experimental and cosmological constraints 

In order to translate the spectra defined by the soft supersymmetry breaking terms in Sections 2.1 
and 2.2, a low energy soft Lagrangian is obtained by solving the renormalization group equations 

■''In fact, the model considered in this reference used the hnear multiplet forrnahsm for the dilaton in which such 
nonperturbative corrections are more easily incorporated into the low energy effective supergravity theory. For a 
description on how to correctly translate from one formulation to the other, see Appendix A of [6]. 



14 



(RGEs) from some initial boundary condition scale to the elcctrowcak scale. Having thus extracted 
the soft supersymmetry breaking parameters at the low-energy scale, the physical mass spectra for 
the superpartncr and Higgs sector can be obtained and various direct collider constraints and indi- 
rect non-collider constraints can be applied to examine the region of phcnomenological viability of 
these models. In this section we describe the tools we use and the various observational constraints 
we impose. 

3.1 Electroweak symmetry breaking and physical masses 

We have used the Fortran code SuSpect [28] to solve the RGEs for the soft supersymmetry breaking 

parameters between the initial high energy scale ^uv and the scale given by the Z-boson mass. 
While the initial scale ^Uuv should itself be treated as a model-dependent parameter, we have 
chosen /Xuv = I^gvt throughout. We use tan/3 and the sign of the supersymmetric /x parameter in 
the superpotential as free parameters, defined at the low-energy (electroweak) scale. 

The magnitude of the /i parameter is determined by imposing electroweak symmetry breaking 
(EWSB) at the scale defined by the geometric mean of the two stop masses (m^^m^^)^^^! which 
minimizes the one-loop scalar effective potential [29, 30]. The one-loop corrected n term is obtained 
from the condition 



where Sm^^ and Sm^j^ represent the one loop tadpole corrections to the running Higgs masses rn^u^ 
and [31, 32, 33]. SuSpect includes the corrections from all the third generation fermion and 
sfermion loops as well as loops from gauge bosons, Higgs bosons, charginos and neutralinos. As the 
superpartncr spectrum depends on the value of the ^ parameter computed in (3.53), the proper 
inclusion of superpartncr thresholds in the KG evolution requires that the procedure be iterated 
until a coherent and stable value for is obtained. Usually, SuSpect requires only two or three 
iterations to have a relative precision of the order of 10~^. We have used three. The soft super- 
symmetry breaking parameters at the weak scale are then passed to the C code micrOMEGAs [34] to 
perform the calculation of physical masses for the superpartners and various indirect constraints, 
to be described in the following section. 

The specific algorithm we employ can be described as follows. The one-loop order SUSY break- 
ing parameters of the heterotic orbifold models obtained in the previous section are entered into 
SuSpect at the high energy scale. Other parameters, including Standard Model fermion masses 
and gauge weak couplings as well as certain EWSB parameters such as tan/3 are input at the low 
energy scale. SuSpect runs the renormalization group evolution of the parameters back and forth 
between the low energy scales such as Mz and the electroweak symmetry breaking scale, and the 
high-energy scale such as the GUT scale. This is the case for the soft SUSY-breaking terms (scalar 






(3.53) 



tan2 (3-1 
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and gaugino masses, bilinear and trilinear couplings and tan P) and fj,. This procedure has to be 
iterated in order to include SUSY threshold effects or radiative corrections due to Higgs and SUSY 
particles. In the first step, these thresholds are only guessed since the spectrum has not been cal- 
culated yet, and the radiative corrections arc not implemented. The thresholds are properly taken 
into account in subsequent iterations. At the elcctroweak scale, the consistency in the calculation 
of the n term is checked using the expression in (3.53). We then call micrOMEGAs, which calculates 
the SUSY spectrum at one loop order using FeynHiggs for the Higgs sector. micrOMEGAs also 
calculates the b .57 branching ratio as well as the relic density of the lightest neutralino and 
the anomalous magnetic moment of the muon. We have checked the consistency between the two 
spectra generated by SuSpect and by micrOMEGAs: the maximum deviation was of the order of 4% 
and was confined to specific regimes such as high tan (3. 

The first condition that we require on a given set of soft supersymmetry breaking masses is 
that an appropriate vacuum state appears - in particular that EWSB occurs. Clearly it is possible 
for the value of /Z^ in (3.53) to become negative for a particular choice of model parameters. In 
such a situation EWSB does not occur. The precise domain of parameter space where this occurs 
is sensitive to the values of the Yukawa couplings - and hence to the choice of quark pole masses 
assumed. This makes the SUSY corrections to the bottom quark mass quite important. We have 
performed the analysis using the heavy fermion masses in SuSpect of 

Mt = 175.0 GeV, Mb = 4.62 GeV, = 1.778 GeV, (3.54) 

and gauge couplings at the Z-mass of 

a|^(M^) = 1/127.938, af^(M;.) = 0.118, = 0.23117. (3.55) 

We also reject all points in the parameter space which give a tachyonic Higgs boson mass (in 
particular, where m\ < 0) or negative sfermion squared masses. 

We next require that the lightest supersymmetric particle (LSP) be neutral, and thus also 
reject all model parameters where the lightest neutralino Xi is not the LSP. In most of the rejected 
cases it is the lightest stau ri or the gluino g which becomes lighter than the Xi- The remaining 
parameter space is further reduced by limits on superpartner and Higgs masses from various collider 
experiments. 

Concerning the mass bounds we use, we take the most recent results combined by the LEP 
Working Group [35, 36]. For the light CP-even neutral Higgs mass (m/i), we assume that a 95% 
confidence level (CL) lower limit on rrih is set at 111.5 GeV. The search for an invisibly decaying 
Higgs boson in hZ production has allowed a 95% CL lower limit on nih to be set at 114.4 GeV, 
assuming a production cross section equal to that in the Standard Model and a 100% branching 
fraction to invisible decays [37]. We believe the value of 113.5 GeV serves as a good reference point, 
and this is the value we use in our analysis. Concerning the chargino limit, we take 103.5 GeV, 
bearing in mind that in some degenerate cases and for light sleptons the limit can go down to 88 
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GeV [38]. For the squark sector the hmit of 97 GeV [39] is implemented. For aU mass bounds 
we should keep in mind that experimental limits are always given in the context of a particular 
SUSY model which is not generally a string motivated one. The bounds we use could possibly be 
weakened in some cases. 

3.2 Indirect constraints on supersymmetric spectra 

Various non-collider observations can be used to further reduce the allowed parameter space of 
these loop-dominated orbifold models. We will focus our attention on the three such processes that 
yield the most stringent constraints: the density of relic neutralino LSPs, the branching ratio for 
decays involving the process h ^ and the measurement of the anomalous magnetic moment of 
the muon. 

3.2.1 The relic neutralino density 

The existence of dark matter is one of the first glimpses of possible physics beyond the Standard 
Model. It is probable that dark matter consists of some stable or extremely long-lived particle left 
over from the hot early universe. For any given extension of the Standard Model containing stable 
or long-lived particles, the present-day relic density of such particles can prove quite constrainting. 
The basic calculation of a relic density from thermal considerations is standard [40] . Here we review 
the steps in such a calculation under the explicit assumption that the lightest neutralino constitutes 
the cold dark matter particle. 

The initial number density of the neutralino is determined because the particle is assumed to 
have been in thermal equilibrium. When the particle begins in thermal equilibrium with its sur- 
roundings, interactions that create neutralinos usually happen as frequently as reverse interactions 
which destroy neutralinos. Once the temperature drops below T ~ m^, most particles no longer 
have sufficient energy to create neutralinos. Now neutralinos can only annihilate, and these anni- 
hilations occur until about the time when the Hubble expansion parameter becomes larger than 
the annihilation rate, H > Fann- When expansion dwarfs annihilation, neutralinos are being sep- 
arated apart from each other too quickly to maintain equilibrium. This happens at the freeze-out 
temperature, usually at Tp ~ for cold dark matter. 

In most neutralino relic density calculations, the only interaction cross sections that need to 
be calculated are annihilations of the type XX ^ ^ where x is the lightest neutralino and X is 
any final state involving only Standard Model particles. However, there are scenarios in which 
other particles in the thermal bath have important effects on the evolution of the neutralino relic 
density. Such a particle annihilates with the neutralino into Standard Model particles and is called a 
coannihilator [41]. To serve as an effective coannihilator, the particle must have direct interactions 
with the neutralino and must be nearly degenerate in mass. Such degeneracy happens in the 
MSSM, for instance, with possible coannihilators being the lightest stau [42], the lightest stop [43], 



17 



the second-to-lightest neutrahno or the Ughtest chargino [44, 45]. When this degeneracy occurs, 
the neutrahno and ah relevant coannihilators form a coupled system. In this section, we will denote 
particles belonging to that coupled system by Xi- Now all interactions involving particles in this 
coupled system come into play, including XiXj ~^ ^7 Xi^ ^ Xj^ ■, s-i^d Xi ~^ Xj^- Here both X 
and Y denote states including Standard Model particles. Decays once again enter the calculation 
because the coannihilators are generally not stable and eventually decay into the lightest neutrahno. 

For the case without coannihilations, evolution of the relic particle number density, n, happens 
in accordance with the single species Boltzmann equation 

^ = -3Hn - {av) [n^ - (n"^)^] , (3.56) 

where n^^ is the equilibrium number density, H is the Hubble parameter at time t, and {av) is the 
thermally averaged annihilation cross section. The number density is modified by Hubble expansion 
and by direct and inverse annihilations of the relic particle. The relic particle is assumed to be 
stable, so relic decay is neglected. In the above expression, we have also assumed T invariance to 
relate annihilation and inverse annihilation processes. 

In the presence of coannihilators, the Boltzmann equation becomes more complicated but can 
be simplified using the stability properties of the relic particle and the coannihilators (using n = 
^i). Application of these simplifications leads to 

dn ^ 

— = - "^Hn - ^ {(TijVij) (niUj - n-'^n^") . (3.57) 

To a very good approximation, one can use the usual single species Boltzmann equation for 
the case of coannihilations if the following replacement is made for the thermally averaged cross 
section: 

For scenarios involving coannihilations, the expression for the thermally averaged cross section 
begins as a six dimensional integral, seven dimensional including the integration over the final state 
angle. This integral has been conveniently put into the form of a one-dimensional integral over the 
total squared center of mass energy [46] . Computer codes exist which numerically perform the center 
of mass momentum integration and the final state angular integration. Some codes include only a 
subset of all possible coannihilation channels [47, 48], while the program micrOMEGAs includes all 
relevant coannihilation channels. Analytical expressions for the cross sections after integration over 
final state angle also exist in the literature [49]. In practice, the Boltzmann equation is generally 
not given a full numerical treatment since accurate approximate solutions exist [40]. 
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Once the effective thermally averaged cross section has been derived as a function of tempera- 
ture, it is used to iteratively determine the freeze-out temperature. In practice, the dimensionless 
inverse freeze-out temperature, xp = m-^/Tp, is calculated using 

xp = ln( 0-038 gmp,m,M\ 

\ y^g^ J 

Here rripL is the Planck mass, g is the total number of degrees of freedom of the x particle (spin, 
color, etc.), 51* is the total number of effective relativistic degrees of freedom at freeze-out, and the 
thermally averaged cross section is evaluated at the freeze-out temperature. For most cold dark 
matter candidates, xp — ^0. 

When the freeze-out temperature has been determined, the total (co) annihilation depletion of 
the neutralino number density can be calculated by integrating the thermally averaged cross section 
from freeze-out to the present temperature (essentially T = 0) . The thermally averaged cross section 
appears in the formula for the relic density in the form of the integral J (xp) = x'^dx: 

V ^pi {g*s/g* ) J {xp) 

Here g^^s is the number of effective relativistic degrees of freedom contributing to the entropy of 
the universe and h is the reduced Hubble parameter, defined by Hq = 100 /i kmsec"^ Mpc~^. The 
above expression is commonly given in the less explicit form: 

".-^ . (3.61) 

gj mpj^J{XF) 

This is the expression that one compares with an experimental determination of the dark matter 
abundance. 

Recent evidence suggests [50] that ~ 0.3 with h^ ~ 0.5. We will take as a conservative 
favored region, 

0.1 < n^^h^ < 0.3. (3.62) 

The lower bound comes from the requirement that Xi should at least form galactic dark matter, 
and the upper bound is a very conservative interpretation of the lower bound on the age of the 
Universe. Recently, WMAP [51] has elegantly confirmed the composition of the Universe to be 
73% dark energy and 27% matter. WMAP determines a total matter density J^mIi^ = 0.135lo:oo9 
and a total baryon density fi^h^ = 0.0224 it 0.0009, from which one can extract a 2a range for 
the density of cold dark matter: r^cDlvlh^ = 0.1126lQgJgJ . We do not use this new result in our 
analysis. The effect of the WMAP measurement has been studied in the case of neutralino dark 
matter from mSUGRA [52] and rSUGRA [53] and also in the case of scalar dark matter from 
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'little Higgs' theories [54]. Let us stress that the requirement of (3.62) should not be treated as 
a constraint, but rather as an indication of the region preferred by cosmological considerations. 
Theoretical assumptions made to extract the a present relic density of neutralino LSP need not 
hold. In fact, the missing non-baryonic matter in the universe may not consist of relic neutralinos 
at all. 



3.2.2 The 6 — ^ 57 constraint 

Another observable where the supersymmetric contribution can be important and measurable is the 
flavor changing decay b ^ sj [55]. In the Standard Model, this process is mediated by virtual isospin 
+1/2 quarks and W^-bosons. In supersymmetric theories, the spectrum allows new contributions 
involving loops of charginos and squarks or top quarks and charged Higgs bosons. As these two 
contributions appear at the same order of perturbation theory, the measurement of the inclusive 
decay B — >■ Xsj is a powerful tool in the search for physics beyond the Standard Model. For our 
analysis, wc will use the results given by the CLEO and BELLE collaborations [56]. 

The Particle Data Group [57] summarizes these results and gives the current limit as: 

BR(b ^ S7) = (3.37 ± 0.37 ± 0.34 ± 0.24;g;f^ ± 0.38) x 10"^, (3.63) 

where the three first errors represent respectively the statistical errors, the systematic error, and 
the estimated error on the model describing the behavior of the quarks in the B-meson decay. 
The fourth uncertainty is the error made by an extrapolation to the entire energy range for the 
photon (cut at 2.1 GeV for the experiment). The last error value is an estimate of the theoretical 
uncertainties. To be as conservative as possible, one could add linearly all the uncertainties of (3.63). 
But we will adopt the procedure taken in the recent benchmark study of Battaglia et al. [58] and 
choose to impose the constraint 

2.33 X lO-'' < BR(b ^ S7) < 4.15 x lO"''. (3.64) 

Let US note that we perform these calculations under the assumptions of minimal flavor violation. 



3.2.3 The muon anomalous magnetic moment 

The relation between the spin s of the muon and its magnetic moment /j, is given classically by 

= 2m where 5^ = 2(1 + a^). Following [59] we introduce the parameter (5^ to quantify the 
difference between theoretical and experimental determinations of a^: 

= (a^ - 11 659 000 x 10"^°) x 10^°. (3.65) 

Recently, the Brookhaven collaboration has given a new measurement of the anomalous mag- 
netic moment of the muon [60] 

i = o^^P = (11 659 202 ± 14 ± 6) X 10"^°, (3.66) 
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where the first error is statistical uncertainty and the second is systematic uncertainty. From this 
the current experimental determination of the parameter 5^ is (5^^p = 203 it 8. 

To understand the implications of this result for supersymmetry one needs to know the Standard 
Model contribution. Unfortunately, the computation of (5™ in the Standard Model is complicated 
by our poor understanding of the hadronic contributions to certain vacuum polarization diagrams. 
The theoretical result depends on whether one uses the e+e~ annihilation cross section or r decay 
data to estimate these contributions [61]: 

^exp _ ^SM ^ ^33 7 ^ -Lj_2) X 10-^° [e+e-] (3.67) 

- = (9-4 ± 10.5) X 10-^° [r-decay], (3.68) 

corresponding to a discrepancy between the Standard Model prediction and experiment of respec- 
tively 3.0 and 0.9 standard deviations. In our discussion we will be less conservative than the 
authors of [59] and consider a 2 standard deviation region about the anomalous moment of the 
muon based on the r decay analysis: 

-11.6 < ^newphysics^ jexp_^SM ^ 3Q 4 ^ (g gg) 

The contribution of SUSY particles to the anomalous moment of the muon mainly comes from 
neutralino-smuon and chargino-sneutrino loop-induced processes. As the chargino-slepton-lepton 
couplings have the same form as the chargino-squark-quark couplings, the parameter space restric- 
tion imposed by the constraint (3.69) will be similar in many respects to the one given by the 
h s'^ process. Thus, the SUSY contribution to (5^ — 2) is large for large values of tan (3 and small 
values of the soft-breaking masses. It is interesting to note that the sign of the SUSY contribution 
is equal to the sign of the /x parameter. With our conventions this implies that positive /x is favored, 
as in the h 57 constraint [55]. 

Let us briefly summarize the above-mentioned constraints that are to be applied below. First, we 
demand that electroweak symmetry be correctly broken, as defined by deriving a positive value for 
Ji^ from (3.53). It should be noted that we use parameters of the Standard Model as defined in (3.54) 
and (3.55). Next, we demand that the mass of the lightest scalar Higgs boson be greater than or 
equal to 113.5 GeV. Additionally, we require that the mass of the lightest chargino be greater than 
or equal to 103.5 GeV. We also require that the lightest slepton mass be larger than 88 GeV and 
the lightest squark mass be greater than 97 GeV, but these two scalar fermion mass limits exclude 
no regions of parameter space in the plots resulting from our analysis. We also indicate regions of 
parameter space that result in the right thermal abundance of relic neutralinos: 0.1 < O^/i^ < 0.3. 
While we technically do not exclude any regions solely for failing to produce the cosmologically 
preferred relic density, we do require that the lightest neutralino be the LSP. In terms of the flavor 
changing decay h S7, we take a conservative approach and require the branching ratio to obey 
2.33 X 10^^ < BR(b 37) < 4.15 x 10~^, given in equation 3.64. Finally, wc follow the recent r 
decay analysis of the muon anomalous magnetic moment and require -11.6 < physics ^ 3Q_4^ 
given above in equation 3.69. 
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4 Results and analysis 



We are now in a position to investigate the way in which low energy observations can distinguish 
different regions of the weakly-coupled heterotic string parameter space. We first recall the stan- 
dard analysis in the context of the minimal supergravity (mSUGRA) parameter space. While the 
unified mSUGRA paradigm is unlikely to find a manifestation in realistic string-based models, it 
nonetheless provides a useful benchmark for calibrating results from one study to the next. 

4.1 Revisiting mSUGRA 

The minimal supergravity paradigm is based on the assumption that the soft supcrsymnictry break- 
ing Lagrangian is determined by only five parameters. These include a universal gaugino mass Miy2, 
a universal scalar mass Mq and a universal trilinear soft parameter Aq, as well as the sign of the 
parameter and tan f3. A quick survey of the soft terms presented in Section 2 indicates that there is 
no point in the parameter space of the string models that we consider that generates such a univer- 
sal outcome. A unified scenario is possible only in the event that the tree-level dilaton domination 
case arises [21], though no explicit realization of this outcome exists in a complete model. Even 
if such a scenario were realized, the soft terms would be constrained to a particular point of the 
mSUGRA parameter space such that 

-^0 = Ml/2 = VSMo. (4.70) 

Despite this lack of theoretical motivation, the extreme simplicity of the mSUGRA approach 
makes it an attractive scenario for phenomenological study [62]. We will here just present the main 
outcomes of our analysis. There are still large areas of parameter space in the {M1/2, Mq) plane 
which are allowed, particularly for large tan/3 where the Higgs mass bound is less constraining. 
The constraint from the process 6 — ^ 57 can be the dominant constraint for low gaugino masses. 
Nevertheless, it is never very severe for positive //, even for high tan/3 (note the change in scale in 
Figure 2 relative to Figure 1). The situation would change in the case of large trilinear A-terms 
because of the possibility of light third generation squarks [43]. However, in this case charge and/or 
color breaking (CCB) problems often arise. 

In both plots of Figure 1 the exclusion regions for the chargino and Higgs masses arc given. These 
exclusion regions implement the experimental limits of 103.5 and 113.5 GeV, respectively, that were 
discussed in Section 3.1. We also provide contours of constant Higgs mass of m/j = 115.5 GeV and 
117 GeV (the dashed lines from left to right) as well as a contour of chargino mass m ± = 500 GeV 
(dotted line). For small and moderate values of tan/3 the requirement that the density of the 
lightest neutralino accounts for the dark matter density is very constraining and sensitive to the 
precise value of tan (3. For moderate and low values of this parameter, one usually finds too much 
relic density in the mSUGRA model. For example the "bulk" region at low scalar and gaugino 
masses is excluded by the Higgs mass limit. Only three regions survive: a narrow band along the 
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Figure 1: Constraints on the mSUGRA parameter space for tan/? — 5 (left) and tan/? = 35 (right). 

Constraints on the (A/1/2, Mq) mSUGRA plane are given for /i > and Aq = 0. The red contours and shaded region 
represent the 0.1 < Q^h^ < 0.3 preferred region. A small region at very low gaugino masses marked "EWSB" is ruled 
out by improper EWSB. Contours of constant Higgs mass of rrih = 115.5 GeV and 117 GeV (the dashed lines from 

left to right) as well as a contour of chargino mass m ± = 500 GeV (dotted line) are given. For a description of the 

^1 

experimental constraints applied, see Section 3. 



r-LSP exclusion zone, a large Higgs pole region at extremely high tan /?, and the "focus point" [63] 
region at high values of Mq (the thin strip along the "no EWSB" region in Figure 2). The region 
next to the r-LSP exclusion zone provides for efficient tix? coannihilation, considerably reducing 
the neutralino density. For high values of tan /3, interesting areas exist where r^^h^ becomes smaller 
due to near-resonant s-channel annihilation through the heavy Higgs states {A or H). Here 
and niH become smaller and their couplings to the b quark and the r lepton increase. In this region 
XiXi ~^ bb, T^T^ dominates and causes significant depletion of the relic density. This can clearly be 
seen in the plot for tan /3 = 50. It should be noted that these results are strongly dependent on the 
treatment of the radiative correction for the bottom mass and the Higgs masses. As stated earlier, 
we have taken = 4.62 GeV. Small changes in these values have dramatic consequences on the 
density. At large Mq and small there is another region (the "focus point" region) where the 

lightest neutralino as well as the lightest chargino acquire relatively small mass and large Higgsino 
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Figure 2: Constraints on the mSUGRA parameter space for tan/3 = 50. Constraints on the (M1/2, Mo) 
mSUGRA plane are given for /i > and Ao = 0. Note the change in scale relative to Figure 1. The red region 
represents the 0.1 < ^l^Yi^ < 0.3 preferred region. The three shaded regions in the lower left are excluded by (from 
left to right) the chargino mass limit, the lightest Higgs boson mass limit and the 6 ^ 37 rate. 



components. The coupling to the bosons W and Z is large enough to make the annihilation process 
into WW or ZZ efficient. 

Our results are in agreement with the previous studies of mSUGRA [62] . The only discrepancies 
appear in the high tan /5 regimes, mostly due to the different treatment of the 6— quark mass and 
associated Yukawa coupling, and the Higgs sector masses. These differences can be understood in 
the light of the work done by Allanach et al. [64]. Indeed, the treatment of the Yukawa coupling of 
the 6— quark (Ab) in SuSpect results in a higher value for A;, than from the other publicly available 
codes (a 4 to 8 percent effect). This discrepancy affects the Higgs masses {niA in particular) via the 
RGEs which are dominated by Af, for high tan /5. For instance, for tan /5 = 50, we can easily find a 
20 to 50 GeV difference for niA, explaining differences in the parameter space allowed by the relic 
density in the A— pole sector. In the remainder of our paper we will restrict ourselves to tan f3 < 35. 
In this range, the discrepancies between codes are smaller than 3 percent. In terms of relic densities, 
the differences between results using the codes DarkSusy [48] and micrOMEGAs [34] come from both 
the 6-mass treatment and the selection of included coannihilation channels: previous versions of 
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DarkSusy did not include any coannihilation channels with sleptons and scalars while micrOMEGAs 
does.^ 

To summarize, in the case of the mSUGRA model, even if large areas of the (M1/2, Md) param- 
eter space are still allowed by present accelerator data, combining the 6 ^ 57 limits, the neutralino 
relic density constraints and the Higgs mass limit selects the r coannihilation region, the "focus 
point" region, or extremely high values of tan /?. These three regions represent the extreme bound- 
aries of mSUGRA. Thus, it is natural go beyond mSUGRA to see how these constraints affect 
orbifiold models. 

4.2 Moduli domination: the anomaly-mediated limit 

The phenomenology of the original (minimal) AMSB model has been extensively studied (for some 
of the earlier work, see [65]). However the phenomenology of the string-based PV-AMSB model, 
defined by the soft terms given in (2.48), has yet to be thoroughly explored. We view this anomaly- 
dominated regime as a particular limit of string models where supersymmetry is broken by the 
F-term of some compactification modulus. Note that this anomaly-mediated limit is one in which 
all tree level soft supersymmetry breaking terms vanish. We thus choose to begin our study of 
string-based heterotic orbifold phenomenology at one loop with this particularly simple regime. 

Like the minimal AMSB scenario, the PV-AMSB model has two free parameters, apart from 
tan/? and the sign of the /i term. One parameter is the overall scale, given by the gravitino mass 
M3/2. The other parameter is the regularization weight p which determines the size of the scalar 
mass terms relative to the gaugino masses and A-terms. In this sense the weight p plays a role 
comparable to the bulk Mq postulated in the minimal AMSB scenario but it is not an ad hoc 
parameter and it appears elsewhere in the spectrum. Note however that the gaugino masses are 
independent of this weight. Thus achieving a sufficiently heavy chargino implies a lower bound on 
the absolute scale as displayed in Figure 3. 

If we look at the branching ratio b ^ we observe than it becomes important in two situations: 
small values of the gravitino masses and/or large values of p. As previously mentioned, small 
gravitino masses give small values of m ± , which leads to a big contribution of the diagram involving 

X^q in the loop for b — >■ 57. Having a large value for p also leads to small squark masses and thus a 
large contribution from the squark exchange diagrams. Moreover, the same limit p ^ 1 also leads 
to light X3,4) X2j which can now also contribute substantially to 6 57 (because of the small /x 
parameter) . 

Concerning the Higgs mass, our findings agree with intuition gained from mSUGRA. One needs 
significant loop contributions from relatively heavy squarks to satisfj- the LEP constraint, especially 
at low tan/3. Thus regions where squarks are light arc ruled out. This includes the region with 
p — > 1 and smaller gravitino mass scales. In Figure 3 wc show the region excluded by our Higgs 

®The upcoming version of DarkSusy will include these coannihilation channels. 
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Figure 3: Constraints on the PV-AMSB parameter space for tan/3 = 5 (left) and tan/3 = 35 (right). 

Constraints on the {M'j,/2,p) plane are given for ^ > 0. The regularization weight p can be as large as p = 1, at 
which point scalar masses vanish at the one loop level (the standard AMSB limit). The areas marked My^ and Mh 
are ruled out by the chargino mass bound and the Higgs mass bound, respectively. The area marked 6 ^ 57 gives 
Br(6 57) > 4.15 X 10"* (a small region where this bound is relevant occurs near p — 1 for tan/3 — 5). The 
horizontal (dotted) contour is a constant chargino mass of m ± = 130 GeV. 



mass constraint in the left plot, as well as a contour of nih = 114 GeV for comparison. In the right 
plot for tan /3 = 35 there is no constraint in this gravitino mass range arising from the Higgs search 
limit and we provide contours of rrih = 117 GeV and 118 GeV. 

In previous AMSB studies it was found that the neutralino thermal relic density is generically 
too small to explain the amount of dark matter [66]. Due to the low ratio of M2/M1 in the AMSB 
scenarios, the wino content of the LSP is quite high. Additionally, coannihilation between the LSP 
and the lightest chargino is also very efficient. Both of these effects combine to make the thermal 
relic density of LSP negligible. Throughout the parameter space of Figure 3 the thermal relic 
density is around ~ 10~^. Thus the anomaly-mediated character of the gaugino sector in this 
model necessitates a non-thermal production mechanism for neutralino LSPs, or another candidate 
for the cold dark matter must be postulated. 
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Figure 4: Constraints on the PV-AMSB parameter space for A/3/2 = 40 TeV. Constraints on tiie (p, tan/3) 
plane are given for /i > 0. Regions excluded by EWSB failure (top), excessive rate for 6 ^ S7 (right) and the Higgs 
mass limit (bottom) are shaded. We have labeled contours of rrih ~ 115.5 GeV and 117 GeV (dashed lines), as well 
as unlabeled contours of rn^± — 112.5 GeV and 113.5 GeV (dotted lines). The solid contour in the lower left is a 
contour of BR(6 37) = 3.53 x 10*. 



We sum up these constraints in Figure 4 where we have shown the allowed region in the {p, tan /?) 
plane for a gravitino mass of = 40 TeV. We clearly see that high tan/3 area is excluded by 
b — > 57, and the low tan f3 area is disfavored because of the Higgs mass limit. In every case, p 
cannot be too large (must be less than ~ 0.85). This excludes the "minimal" AMSB scenario 
{p = 1) up to the two-loop corrections to soft terms which have not been fully calculated in 
these string-based models. Most of the interior of the parameter space in Figure 4 has sufficiently 
heavy Higgs masses, as shown by the dashed contours of nih = 115.5 GeV and 117 GeV. The 
supersymmetric contributions to the b ^ rate in the interior are small as well. The solid 
contour is BR(6 — > 57) = 3.53 x 10"'^. The chargino masses are very light throughout and we 
have provided contours of m ± = 112.5 GeV (right-most dotted line) and 113.5 GeV (left-most 
dotted line). We note that much of this parameter space predicts a relatively small supersymmetric 
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contribution to the muon anomalous magnetic moment due to relatively heavy scalar masses {i.e. 
S^^^^ ~ 0). The upper shaded region is ruled out by improper electroweak symmetry breaking 

< 0). 

In summary, the Pauli-Villars rcgularization weight p generates a whole class of anomaly- 
dominated models at the one-loop level. While this new degree of freedom solves the tachyonic 
slepton problem, the minuscule relic density still does not allow the neutralino to be a thermal 
dark matter candidate. Additionally, the constraints coming from accelerator physics exclude high 
values of p and low values of M3/2 and predict small contributions to the muon anomalous magnetic 
moment. 

4.3 The general moduli dominated case 

We will now relax the assumption that the various Kahler moduli are stabilized at self-dual points, 
while maintaining the assumption that the modular weights of all matter fields are equal and given 
by = —1- When the vacuum value (Re t) is not fixed at (Re t) = 1 we must use the soft 
terms in (2.47). Note that the value of the Green-Schwarz coefficient 5qs becomes relevant for the 
determination of gaugino masses in this case. For the time being we will take Sqs = in order to 
keep the parameter space to manageable size. 

For the study of the moduli dominated case we have looked at the behavior of the model as 
a function of the vacuum expectation value (vev) of the real part of the (universal) modulus field 
t = T\0=o, for different values of p and tan/3. To be more precise, we have taken two values of p, 
and 0.95, and two values of tan/3, 5 and 35. We have chosen to use the value p = 0.95 because at 
this value the gaugino masses and scalar masses have roughly the same magnitude. It is also very 
near the pure AMSB limit. We first stress the common features of the moduli-dominated case, and 
extract the corresponding phenomenological consequences, before looking at specific points in the 
parameter space. 

We choose to investigate the region 0.1 < (Rc t) < 5. Since the vacuum value of this modulus is 
related to the radius of the compact dimensions, negative values have no physical meaning. Let us 
note that the combination [t + t)G2 {t, t), which appears in the gaugino mass when the equation of 
motion (1.11) is substituted into (2.47), is invariant (up to a sign) under the duality transformation 
t ^ 1/t as shown in Figure 5. However, since this term competes with the anomaly-mediated 
contribution to gaugino masses at the one loop level, gaugino masses will not display this duality 
symmetry. Any value (Re t) 1 represents a spontaneous breaking of modular invariance and the 
phenomenology of the theory for (Re t) < 1 will be different from the case (Re t) > 1. 

Looking again at the gaugino masses in (2.47) it is clear that when 5qs = the gaugino masses 
arc proportional to beta-function coefficients. Thus we expect the phenomenology of this scenario 
to be very similar to that of the PV-AMSB scenario of the previous section. For example, we 
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continue to have 



Ml 



gi{^GVT)h h = E=6s (4.71) 



GUT 52(MGUt) b2 62 5 

independent of the value of (Re t) or p. At the low scale this impUes M2 -C Mi and the lightest 
neutralino Xi and the lightest chargino xf are in a nearly complete Wino state. This means that for 

all values of p and (Re t) we expect m o ~ m ± ~ M2, and the relic density of the LSP neutralino 

is almost entirely depleted by (co)annihilation channels involving Xi cind xf- 

Unlike the case of Section 4.2, however, there is now a value of (Re t) such that the contribution 
from the Eisenstein function exactly counterbalances the contribution from the superconformal 
anomaly giving us vanishing gaugino mass soft breaking terms M^ to this order. For the choice of 
phase conventions in (1.7) and (1.8) this occurs when^ 

{t + t)G2 {t,t) = l C{t) = ^ Ret = 0.523. (4.72) 

Not surprisingly, then, we find that much of the parameter space in the vacinity of (Re t) = 0.5 is 
ruled out by the chargino mass bounds for any reasonable choice of scale M3/2 in Figures 6 and 7. 

In the extreme cases where (Re t) ^ 0.1 and (Re t) ^ 5 the absolute value of the combination 
(t + t)G2 {t,t} becomes large and dominates over the anomaly contribution in the gaugino masses. 
Here gauginos will typically be similar in size to the scalar masses at the high scale, particularly as 
the value of p is increased towards its limiting value p = 1 (Figure 7) . The chargino mass bound is 

'^In [6] the same behavior was noted. In that reference, however, the opposite sign convention on (1.7) was used, 
leading to vanishing gaugino masses when (t + t)G2 {t,t) = —1 which occurs at the dual point Re t ~ 2. 



29 



tanp = 5 tanp = 35 




1.0 2.0 3.0 4.0 5.0 1.0 2.0 3.0 4.0 5.0 

< Re T > < Re T > 



Figure 6: Constraints on the moduli-dominated parameter space for tan/? = 5 (left) and tan/3 = 35 
(right) with p — 0. Constraints on the (M3/2, (Re t)) plane are given for Scs = and fi > 0. The region in (Re t) 
shown here is 0.2 < (Re t) < 1. The areas marked and Mh are ruled out by the chargino mass bound and the 

Higgs mass bound, respectively. Chargino masses of m ± = 200 GeV (dotted line) are provided in both plots, while 

^1 

Higgs contours of irih = 115 GeV and 118 GeV (dashed lines) are given where applicable. 



easily satisfied here and the increased contribution from M3 in the RG evolution of the scalar masses 
makes the Higgs mass constraint easier to satisfy, particularly for large tan (3. In Figure 6 we show 
the Higgs and chargino exclusion regions in the tan/5 = 5 panel and the contours mh = 115 GeV 
(dashed line) and m ± = 200 GeV (dotted line). The Higgs mass is not constraining at tan/3 = 35 

Xi 

and we provide contours of nih = 115 GeV and 118 GeV (dashed lines) and again m ± = 200 GeV 
(dotted line). 

While the mass contours and exclusion regions in Figure 6 show the underlying modular sym- 
metry t ^ 1/t in a very broad sense, it is clear that these contours do not obey this symmetry 
in the strict sense. We have indicated the self-dual point (Re i) = 1 by a heavy vertical line. 
The exclusion regions do not show the expected duality symmetry about this point. As already 
mentioned above, this is due in large part to the competition between contributions to the gaugino 



30 



tan[3 = 5 tanp = 35 




< Re T > < Re T > 



Figure 7: Constraints on the moduli-dominated parameter space for tan/? = 5 (left) and tan/3 = 35 

(right) with p — 0.95. Constraints on the (M3/2, (Re t)) plane are given for Sgs ~ and /i > 0. The shaded region 

in the upper right of the tan /3 = 35 plot is excluded by CCB constraints, while the snail shaded triangle in the 

lower left is excluded by improper EWSB. Chargino mass contours of m ± = 200 GeV (dotted line) and Higgs mass 

^1 

contours of mu = 115 GeV (labeled dashed line) are given in both plots. 



masses from the moduli sector and contributions from the superconformal anomaly. The relative 
sign between these contributions depends on which side of the self-dual point the (overall) modulus 
is stabilized. Having chosen the conventions of (1.7), (1.19) and (1.20) this implies a zero in the 
gaugino mass formula (1.30) on a definite side of the self-dual point (with the conventions here, the 
smaller value side). The conventions can be changed to make this zero occur on the other side, but 
only at the expense of changing the convention for the sign of the ^ parameter. As we keep this sign 
always fixed at positive values, the relative sign between the gaugino masses - in particular, the 
gluino mass M3 - and the ^ term is the key variable in distinguishing the two sides of the self-dual 
point. This is reflected in the fact that at tan (3 = 35 in Figure 6 the b ^ S'j constraint applies to 
only one side of the figure. Other observables, such as the muon anomalous magnetic moment, will 
be sensitive to this relative sign as well, making any spontaneous breaking of modular invariance a 
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true observable, at least in principle. 

As we approach the minimal AMSB limit of p — ^ 1 the relevant constraints change as scalar 
masses diminish. This is displayed in Figure 7. The Higgs mass limit is important over much of 
the parameter space, requiring gravitino mass scales well in excess of 50 TeV for (Re t) ~ 1 in 
order to compensate for the rapidly diminishing stop masses. Note that Figures 6 and 7 correctly 
reproduce our previous results in the limit of (Re t) — ^ 1. The light squarks in the p ~ 1 regime 
tend to contribute too much to the 6 — 57 process - a process that provided no such constraint in 
this region of (Re t) in Figure 6 due to the relatively large scalar masses in those cases. 

Nearly the entire region of the parameter space for tan/3 = 35 shown in Figure 7 is ruled 
out by this constraint. The remaining slender allowed region near (Re t) ~ 0.25 allowed by the 
experimental constraints is the result of a relatively light charged Higgs that helps to cancel the 
contributions from the light chargino. Here the scalars and the gauginos are light - in particular, 
the gluino has a vanishing soft mass at the high scale near (Re t) ~ 0.3 when Sqs = (see Figure 8 
below). This leads to a very small /x term value in this narrow range, and consequently a small 
charged Higgs mass. For example, as we move from (Re t) ~ 0.1 to (Re t) = 0.3 along the line 
M3/2 = 15 TeV in the right panel of Figure 6 the value of n at the electroweak scale drops from 
1700 to 1100 GeV due to the diminishing gluino mass. But in the right panel of Figure 7 n goes 
from 1700 to 380 GeV over the same range, leading to a charged Higgs mass of approximately the 
same value. Scalar masses are around 200 GeV at this point. 

In the right panel of Figure 7 with p ^ 1 (light squarks) and large tan (3 (large couplings in the 
Higgs sector) there also appears a zone excluded by too large of a super symmetric contribution to 
the anomalous magnetic moment of the muon (a black region labeled "/x" ) due to the presence of 
very light squarks and gauginos. This region is already excluded, however, by multiple observables 
such as the 6 ^ 57 rate, the chargino mass and the Higgs mass bound. The shaded region in the 
upper right of the tan/? = 35 plot in Figure 7 is disfavored by the presence of CCB minima deeper 
than the preferred electroweak vacuum. The presence of such a region is due to the very large 
trilinear couplings for these values of (Re t) which are larger than the scalar masses when p 2± 1. 

In conclusion we see that in this simple moduli-dominated scenario the phenomenology is similar 
to that of anomaly mediation. The closer one approaches the minimal AMSB limit of p = 1 the 
more problematic the scenario becomes, particularly at large tan (3. Other points in the PV-AMSB 
class have some areas of phenomenological viability, though the chargino and Higgs mass bounds 
can be quite constraining, particularly near the self-dual point (Re t) = 1. All of these cases will 
be unable to explain the cold dark matter content of the universe, at least in terms of thermal relic 
neutralinos. 

Our study of the PV-AMSB scenario of Section 4.2 had one effective parameter, the regulariza- 
tion weight apart from the overall scale given by the gravitino mass M3/2. This regularization 

*For a discussion of these points, in the context of distinguishing anomaly mediation from gauge mediation, see 
the work of G. Kribs in [65]. 
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weight controls the size of the scalar mass relative to the fundamental scale M3/2 while the gaugino 
masses were fixed by their beta-function coefficients. In the present section the parameter space 
has been expanded to allow for (Re t) ^1, allowing the gaugino masses to be varied relative to the 
scalars and gravitino mass M3/2 in an independent way. Small, or even vanishing, gaugino masses 
are now possible, as is the possibility of relatively large gaugino masses as one moves far from the 
self-dual point for the compactification moduli. 

Despite this new degree of freedom the parameter space is still quite constrained. In particular, 
obtaining the correct relic density for the neutralino LSP is impossible in this framework, because of 
its extreme wino-like nature. But, as we will see in the next section, the possibility of non- vanishing 
Green-Schwarz counterterm can open a new region of parameter space completely in accord with 
present experimental data. 

4.4 The influence of the Green-Schwarz counterterm 

In this section we wish to introduce the possibility that 5^^ 7^ 0, as would typically be the case in 
realistic heterotic orbifold constructions. As we will see, this parameter allows one to interpolate 
between regions with a phenomenology similar to anomaly mediated models (specifically the PV- 
AMSB scenario of equation 2.48) and regions where the phenomenology resembles that of the 
minimal supergravity paradigm. 

In the conventions defined by the string threshold correction (1.28) the value of (5gs is a negative 
integer between and -90. For any given orbifold construction the precise value of this coefficient can 
be worked out, and indeed such an exercise has been performed for realistic Z3 models already [67]. 
In the figures that follow we will treat this variable as a continuous parameter. At values of the 
Kahler modulus far from its self-dual point, where the Eisenstein function is not negligible, the 
term in (2.47) proportional to 8qs can give rise to substantial (universal) contributions to the 
gaugino masses. Thus points in parameter space with large values of 5qs and values of (Re t) 
far from unity will likely have a pattern of gaugino masses similar to mSUGRA - and hence arc 
more likely to provide a viable dark matter candidate than cases with ~ 0. Furthermore, once 
\5qq\ ~ vT&Tr^ ^ C'(IO), the gaugino masses and scalar masses will be of approximately the same 
size. This is very different from the hierarchical situation of the previous sections. 

Because of the difference of sign between the beta functions ha of the Standard Model gauge 
groups, M3 has a particular behavior that differs from the other gaugino masses. This in turn 
drives the phenomenology of the model, particularly through the gluino's RG effects on the rest of 
the superpartner spectrum. With a non-vanishing Green-Schwarz coefficient the various gaugino 
masses will no longer have zeros at the same value of Re t. In particular, for certain combinations 
of (Re t) and 6qs it is possible for the gluino mass to vanish at the boundary condition scale while 
the other gauginos have non-vanishing masses. These combinations are shown in Figure 8. For 
parameter choices near these contours we might expect the phenomenologically unacceptable result 
of a gluino LSP. Even when the gluino is not the LSP, its small value relative to other gauginos may 
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Figure 8: Contours of relative running gaugino masses M3/M2 in the {(Ret), 5c,s) plane. These soft 
masses are at the initial (GUT) scale. The heavy (dark) contour is the limit of vanishing gluino mass (there is 
another such contour in the upper left corner on the other side of the self-dual point). For (Re t) > 1 we also give 
contours of IM3/M2I = 0.33 (dashed) and 0.75 (solid). 



provide too little radiative correction to squark masses and the Higgs mass bound may be difficult 
to satisfy. On the other hand, a relatively light gluino is likely to reduce the amount of fine-tuning 
required to obtain Mz = 91.2 GeV [68]. 

In Figures 9, 10 and 11 we look at this model in the {Sqs, -^3/2) plane for p = and (Re t) = 0.5, 
(Re t) = 1.23 and (Re t) = 2, respectively. Beginning with the case (Re t) = 0.5 in Figure 9 we see 
that the limit as Sqs reproduces the cases studied in the previous section. In particular, such 
models require very large M3/2 in these regions to satisfy the chargino and Higgs mass constraints. 
This behavior is seen in all three of the (Re t) presented here. A maximum value of \Sgs\ can in 
general be obtained by requiring that lightest neutralino be heavier than the lightest stau. For any 
values of tan /3 the maximum value of | (5gs I is determined by the requirement of keeping the lightest 
neutralino as the LSP. The quantity I^qI^I is smaller for higher tan/3 because the stau is in this 
case lighter. 
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Figure 9: Constraints on the moduli-dominated parameter space for tan/3 = 5 (left) and tan/3 = 35 
(right) with p — and (Re t) — 0.5. Constraints on the Sgs) plane are given for /i > 0. Shaded regions to 

the left are ruled out by improper EWSB or a stau LSP. The dark region in the lower right of the tan/3 — 35 plot 
labeled "/i" has too much supersymmetric contribution to (y^ — 2). The strip labeled " f2" has the cosmologically 
preferred neutralino relic density. Dashed contours of — 118 GeV (left) and irih = 122 GeV (right) are also given. 



In between these disallowed regions the increasing values of the gaugino masses for fixed scalar 
mass, as the Green-Schwarz coefficient increases in absolute value, results in smaller contributions 
to the b s-f branching ratio and a smaller (in absolute value) anomalous magnetic moment of the 
muon. In Figure 9 we have chosen a value of (Re t) near the value where gaugino masses vanish 
when 6qs — > 0. As is evident from the figure, gaugino masses are indeed falling as this limit is 
reached. Note that in this case there is no point in the plane for which the gluino is the LSP. 

Concerning the relic density, it is negligible for small values of the coefficient for the Green- 
Schwarz counterterm, as we have seen previously: the neutralino is extremely Wino-like in nature 
and is nearly degenerate with the lightest chargino. Almost all of the relic neutralinos are depleted 
through coannihilation channels with the charginos. When we increase I^gsI) however, we increase 
the gaugino masses through the increased importance of the G{t,t) terms in (2.47). This universal 
contribution can come to compete with the (nonuniversal) anomaly-mediated term, driving the 
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Figure 10: Constraints on the moduli-dominated parameter space for tan/3 — 5 (left) and tan/3 = 35 
(right) with p = and (Re t) — 1.23. Constraints on the (Af3/2, Sas) plane are given for jj, > 0. The lightly shaded 
region in the center gives a gluino LSP. The thin strip with acceptable thermal neutralino relic density is labeled 
"r2." Note the change in scale from Figure 9. 



gaugino mass terms Mj=i^2 far away from the previous "degenerate" situation. The gaugino mass 
sector then begins to look similar to that of mSUGRA and a more bino-like LSP can develop, 
leading to an increased relic density. Just as in the mSUGRA case there continues to be a region 
of stau coannihilation near the excluded stau-LSP area. And - again, just as in mSUGRA - much 
of the "bulk" area is ruled out by the chargino and Higgs mass bounds, or the b ^ s'j constraint. 

In Figure 10 we investigate the same parameter space for the case where (Re t) = 1.23. When 
multiple gaugino condensates are utilized to stabilize the dilaton, with the tree level Kahler po- 
tential given by K{S, S) = — ln(5 + S), then a modular invariant treatment of the resulting non- 
perturbative potential for the moduli fields will lead to the conclusion that the auxiliary field for 
the dilaton must vanish in the vacuum. The remaining potential for the Kahler moduli leads to 
their stabilization at this value [69]. Specific points in the space shown in Figure 10 were singled out 
for a more detailed study of their collider signatures in a recent set of string-inspired "benchmark" 
models [70]. Here a more complete survey of the parameter space is possible. 
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Figure 11: Constraints on the moduli-dominated parameter space for tan/3 — 5 (left) and tan/3 = 35 
(right) with p = and (Re t) — 2.0. Constraints on the {l\I-i/2, <5gs) plane are given for /i > 0. The dark shaded 
regions on the left have a stau LSP. The tan (3 = 35 plot also has a region with a gluino LSP. For tan /3 = 5 the region 
labeled "fi" has the cosmologically preferred relic density of neutralinos. No such region exists for the higher tan P 
plot. In that case the exclusion contours are due to (from bottom right to upper left) CCB vacua, the chargino mass, 
too large SUSY contributions to (gjj — 2), the Higgs mass limit and too large a 6 — > 57 rate. 



In accordance with Figure 8 we note a region centered about Jqs = ~15 where the gluino is 
the LSP. Interestingly, this is precisely the region of values favored by semi-realistic Z3 orbifold 
constructions [67]. To the right, as the absolute value of the Green-Schwarz counterterm diminishes, 
the chargino mass bound is eventually reached. The value (Re t) = 1.23 is sufficiently close to the 
self-dual value that the theory near 6qs = is anomaly mediated-like in character: very large values 
of the gravitino mass are generally required to achieve a sufficiently heavy lightest chargino and 
lightest Higgs (note the difference in scale between Figures 9 and 10). For tan/3 = 35 the Higgs 
mass constraint is less important, being only relevant for low gravitino masses near the gluino LSP 
region. In addition, the LSP is predominantly wino-like on the right side of the gluino LSP region, 
and the relic density of these neutralinos is negligible. 

As we move to larger absolute values of 5qs we interpolate between this anomaly mediated 
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Figure 12: Constraints on the moduli-dominated parameter space for tan/3 — 5 (left) and tan/9 = 35 
(right) with p = and A/3/2 = 10 TeV. Constraints on the {(Ret) , Scs) plane are given for ^ > 0. The upper 
light shaded region has improper EWSB, while the lower light shaded regions on either side of (Re t) = 1 have a stau 
LSP. All collider constraints have been combined in the dark shaded region to make the preferred relic density strip 
easier to distinguish. 



region and an mSUGRA-like region. For certain critical values of this parameter the wino content 
of the LSP is sufficiently reduced to result in the correct annihilation efficiency to produce the 
cosmologically preferred relic abundance. This change in the character of the LSP is the result 
of the increasing importance of moduli contributions to the gaugino masses (the first term in the 
gaugino mass parameter in equation 2.47) versus the anomaly mediated contributions (the last 
term in that expression). The ability to achieve the right relic density, independent of the scalar 
mass values, is a general property of models with this sort of non-degeneracy among gaugino mass 
parameters [44, 45]. 

We investigate the case (Re t) = 2.0 in Figure 11. The behavior at small tan/3 is not dissimilar 
from the nominally "dual" case of (Re t) = 0.5 in Figure 9. The stau LSP constraint now arrives 
at much smaller values of the Green-Schwarz coefficient. Again, the region with the promising relic 
density of neutralinos is largely constrained by the Higgs mass limit. 
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The parameter space for the higher tan/3 regime is even more tightly constrained (at least at 
these relatively low values of the gravitino mass). We again see a region of gluio LSP - a narrower 
region at lower absolute values of ^gs consistent with Figure 8. At the bottom of the plot the values 
are ruled out by the presence of CCB vacua while the region to the left has a stau LSP. The light 
charginos for this value of the gravitino mass lead to large contributions to the muon anomalous 
magnetic moment (the black region labeled "//" ) . Finally, all but the uppermost corner of the space 
is ruled out by the Higgs mass limit and the upper bound on BR(6 — 57) . 

We combine all of these observational and theoretical constraints to explore the relation between 
5gs and (Re t) in Figure 12. We plot the allowed parameter space in the {6gs, (Re t)) plane in a 
similar fashion to Figure 8 for the specific gravitino mass of M3/2 = 10 TeV. On either side of 
the self-dual point at large negative values of Sqs the stau is the LSP. Near these regions there is a 
strip of preferred neutralino relic density. The area at small values of \6qs\ and near (Re t) = 1 is 
ruled out by improper EWSB. Here again there is a small strip of preferred relic density, though 
this is within the region ruled out by the Higgs mass, chargino mass or 6 — 57 constraint. In 
the interior there continues to be viable parameter space, including points with acceptable relic 
neutralino densities. 

4.5 The generalized dilaton domination case 

If we now consider the other extreme case {9 = n/2) the general features are somewhat different. 
Referring to the soft terms of equation 2.52, we can clearly see that the soft scalar mass terms will 
be roughly dominated by the gravitino mass scale (up to radiative corrections). We anticipate, 
therefore, that the relevant gravitino mass scale will be approximately one order of magnitude 
lower than the one necessary in the moduli-dominated case because of the absence of the loop 
suppression factors on the scalar masses. The gaugino mass soft supersymmetry breaking terms 
will be determined by the dilaton auxiliary field vev ( ) . 



In the case where nonpcrturbative corrections to the dilaton Kahler potential are imagined, 
the minimum of the combined (modular invariant) dilaton/Kahler moduli potential now occurs at 
i^P^'^ 7^ and compactification moduli stabilized at self-dual points where ^P'^'^ = 0. As described 
in Section 2.2, the requirement of vanishing vacuum energy naturally leads to a suppression of 
this auxiliary field vev relative to the auxiliary field of the supergravity multiplet. Thus the two 
contributions to the gaugino mass in equation 2.52 each involve a loop suppression factor (unlike the 
scalar masses) and the dilaton and conformal anomaly contributions are comparable. Indeed, the 
beta-function coefficients of the Standard Model gauge groups are of the order of 10~^. Moreover, 
the term will be driven by b+, the largest beta-function coefficient among the condensing gauge 
groups of the hidden sector, from (2.51) and (2.49). In fact, if we look in more detail at the 
expression for F^ in equation 2.51, it is apparent that for not too large values of 6+ we can 
consider that F^ has a linear evolution as a function of that parameter. Increasing b-^- has a 
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tan[3 = 5 tan(3 = 35 




Figure 13: Constraints on the dilaton-dominated parameter space for tan /3 = 5 (left) and tan /3 = 35 
(right) with . Constraints on the (M3/2, &+) plane are given for fi > 0. 



direct consequence on the values of the gaugino breaking terms, while increasing M3/2 has direct 
consequences on the general size of all soft breaking terms - in particular soft scalar masses. We 
can consider the effect of the parameter as a "fine structure" on top of the gross feature of a 
hierarchy between scalar s and gauginos in this regime. 

The general features of this situation are shown in Figure 13. When the suppression of the 
gaugino masses is large from small values of 64. (small a^p), the gaugino sector looks increasingly 
like that of anomaly-mediation. This is true of any value of tan /3. As this value is increased the wino 
content of the LSP diminishes and we approach the overwhelmingly bino-like LSP of mSUGRA. 
However, the value of 6+ cannot increase without limit. As we are considering the weakly coupled 
Eg X £"8 heterotic string, we cannot imagine a condensing group larger than the complete Eg of 
the hidden sector, for which = OO/lGvr^ = 0.57. It is more likely that the same Wilson-line 
mechanism that breaks the observable sector gauge group to the product of sub-groups at the 
compactification scale also breaks the hidden sector group as well. Thus a more typical range for 
this parameter might be b-^- < IS/lGvr^ = 0.095, where this value corresponds to a pure Yang-Mills 
sector of SU(5). 
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For this range of condensing group beta-function, only small values of tan/? are consistent with 
the stringent demands of electroweak symmetry breaking. This fact was observed in this model some 
time ago [71], but it is here displayed explicitly in Figure 13. The large regions marked "EWSB" 
in this figure have /x^ < 0. This is a manifestation of the "focus point," or hyperbolic nature of 
this model [63]: the scalar masses are generally much larger than gaugino masses throughout the 
parameter space. The light gluinos imply that third generation squarks must have large masses 
at the boundary condition scale to ensure a sufficiently heavy Higgs mass. Thus M3/2 ^ 4 TeV is 
needed in the tan /? = 5 case and M3/2 ~ 1-3 TeV is needed in the tan P = 35 case. 

Note also that the very large (and nearly universal) scalar masses will imply heavy sleptons at 
the electroweak scale. Therefore both direct neutralino annihilation through slepton exchange, as 
well as Ti/Xi coannihilation, will be ineffective in reducing the relic neutralino density in the early 
universe. Nevertheless, acceptable relic densities are possible in this model, without the need for 
resonant annihilation through a heavy Higgs eigenstate. The cosmologically preferred region is the 
shaded strip labeled "O" in Figure 13. Between this strip and the region ruled out by improper 
EWSB the relic density is below the preferred amount. In this region the LSP has a sufficiently 
large wino content to annihilate efficiently. Below the marked strip the LSP is bino-like and the relic 
density is too large. We note that the region with the correct relic density can continue indefinitely 
to large scalar masses by maintaining the correct wino-content of the LSP. 

Conclusion and Perspectives 

Even within the restricted context of the weakly-coupled heterotic string compactified on an orb- 
ifold, the phenomenology of such supergravity effective theories is far richer and varied than that 
of the standard "minimal" supergravity approach. While our current indirect knowledge on the 
nature of supersymmetry already constrains these models significantly, there still exist large re- 
gions of parameter space for all the cases studied here - even for low values of tan /3. We anticipate 
that future measurements of, or limits on, superpartner masses and kinematic distributions will 
constrain the parameters of these models further. So too will future cosmological or astrophysical 
measurements. In fact, these sorts of measurements form a useful complementarity that will be 
crucial in unraveling the nature of supersymmetry breaking and transmission in such models in the 
manner begun here. Future studies of string-based models should focus on extending this initial 
survey to true collider signatures in both hadron and lepton machines, as well as computing event 
rates for astrophysical processes. 

Our analysis in this paper was somewhat cursory - we consider the constraints that arise 
principally from the chargino and Higgs searches at LEP as well as the measurement of the branching 
ratio for 6 57 in the minimal flavor violation scenario. We have indicated, where applicable, 
the region that seems to be favored by cosmological observations of dark matter relic densities. 
Interpreting this observation as a constraint on a model implies the assumption that the dark 
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matter is the result of thermal processes in the early universe and is composed entirely of relic 
neutralino LSPs. In string models many other dark matter candidates are known to exist: axions, 
superheavy exotic string states, hidden sector gauge and matter condensates, etc., so care should be 
exercised in applying the cosmological bounds. Wc have therefore kept the range we show within 
the wider region 0.1 < ri^h^ < 0.3 rather than the far more restrictive result from the recent 
WMAP experiment. We have also used only the most conservative possible interpretation of the 
measurement of the muon anomalous magnetic moment as to not prematurely prejudice one model 
versus another. 

Despite the uncertainties that must be borne in mind in such an analysis arising from the 
computational tools employed {e.g. the choice of quark pole masses utilized, uncertainties in the 
determination of p, at large tan (3 etc.) some general statements can certainly be made. Among these 
is the observation that explicit models of moduli stabilization utilizing field-theoretic methods tend 
to predict that at least some sector of the theory obtains soft supersymmetry breaking terms only 
at the one-loop level. This is especially true of the gaugino sector. Furthermore, this suppression 
leads to the general requirement that the gravitino be relatively heavy in these cases - perhaps 
alleviating the cosmological problems so often associated with models that contain gravitinos. 

Parameters that are related to the orbifold can have a large impact on these broad features. 
This is a welcome result - implying that experimental data can indeed probe the nature of the 
underlying theory within a class of models. For example, in cases where supersymmetry breaking 
is transmitted predominantly by the moduli associated with compactification the relative sign of 
terms proportional to their auxiliary field and that of the conformal anomaly can in principle be 
measured. This parameter, in turn, is related to the spontaneous breakdown of modular invariance 
that may occur in such models. 

The phenomenology, even in the simple initial approach taken here, is neither that of minimal 
supergravity nor minimal anomaly mediation but is often a hybrid of the two. This makes a detailed 
study of collider signatures all the more urgent as search strategies are often based on one or the 
other of these paradigms. Given that the era of experimental supersymmetry may well be at hand 
- and that weakly coupled hctcrotic strings still represent the best motivated string-based approach 
to understanding supersymmetry and its breaking in the low-energy world - the time is right for 
this exciting undertaking. 
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